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Abstract 

We develop a new operator quantization scheme for gauge theories in which 
the dynamics of the ghost sector is described by a = 2 supersymmetry. In 
this scheme no gauge condition is imposed on the gauge fields. The corre- 
sponding path integral is explicitly Lorentz invariant and, in contrast to the 
BRST-BFV path integral in the Lorentz gauge, it is free of the Gribov ambi- 
guity, i.e., it is also valid in the non-perturbative domain. The formalism can 
therefore be used to study the non-perturbative properties of gauge theories 
in the infra-red region (gluon confinement). 



Introduction. 



Systems with first class constraints, and in particular relativistic gauge theories, 
posses unphysical degrees of freedom which have to be eliminated in the quantiza- 
tion process It is, however, important to note that in general the process of 
gauge fixing and quantization does not commute as was exemplified in , and that 
the unphysical degrees of freedom have to be eliminated after quantization. The 
conventional procedure of the path integral quantization of gauge theories relies 
on the elimination of unphysical (gauge) degrees of freedom on the classical level 
via gauge fixing and, therefore, should be modified ^. 

In relativistic gauge theories the situation is slightly more complicated since one 
wants to build a Lorentz invariant theory. To this end one has to use Lorentz 
invariant gauges, e.g., the Lorentz gauge. This implies, however, that unphysical 
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bosonic degrees of freedom are dynamically activated and have to be suppressed. A 
general procedure for doing this, based on the BRST symmetry, was given in 0. 

A problem common to all quantization procedures based on imposing a gauge 
condition on the gauge fields, is the Gribov ambiguity |^. This states that it is 
impossible (assuming certain boundary conditions at spatial infinity) to find a gauge 
condition that would eliminate all gauge equivalent configurations ^ . This implies 
that even after gauge fixing there is still a set of residual gauge transformations under 
which the action will be invariant. Although these residual transformations do not 
lead to a further reduction in the number of physical degrees of freedom, they are 
important in determining the topological structure of physical configuration (phase) 
space and therefore they play a crucial role in the quantization process and the 
corresponding functional representation For normal Lorentz invariant gauges, 
such as the Lorentz gauge, this set of residual transformations is extremely complex 
and the task of incorporating them in the quantization process seems hopeless. 

There is good reason to believe that the Gribov ambiguity is irrelevant in the 
perturbative (asymptotic free) domain 0. However, there is evidence that it plays 
an important role in the infra-red, non-perturbative behavior of the gluon propogator 



From the previous introductory remarks it should be clear that one needs to 
develop a quantization procedure which eliminates the unphysical degrees of free- 
dom while it (a) avoids imposing a gauge condition on the gauge fields, i.e., the 
Gribov ambiguity and (b) leads to a Lorentz covariant theory. The basic ingredients 
for such a scheme were given in [12|. The idea of is to extend the theory in 



a supersymmetric and gauge invariant way by introducing bosonic and fermionic 
ghosts transforming in a gauge multiplet. The supersymmetry insures that the con- 
tributions of the boson and fermion ghosts cancel in the partition function, while 
the gauge invariance makes it possible to avoid imposing a gauge condition on the 
gauge fields, but rather to eliminate the unphysical degrees of freedom by imposing 
a gauge condition on the bosonic ghost fields. Finally the whole procedure respects 
the Lorentz invariance of the theory. This program was, however, performed within 
the functional integral setting and, given the ambiguities that arise in the functional 
integral formalism, it is highly desirable to develop this program on the operator 
level and derive the corresponding functional integral representation from there. 

Our aim with the present paper is to develop this operator quantization scheme. 
The steps we follow in doing this are essentially the same as outlined above. We show 
how any quantum mechanical system (regardless whether it has a gauge symmetry 
or not) can be extended by adding N = 2 supersymmetric ghosts. In the extended 
theory physical states are identified as those invariant under SUSY transformations. 
Matrix elements of any system operator calculated in the physical subspace coincide 
with those of the original system. 

In the case of gauge theories this extension is done by putting the ghosts in a 
gauge multiplet. This implies (a) that we not only modify the Hamiltonian, but 
also the constraints and (b) that the extension respects the gauge symmetry, i.e. 
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commutators of the extended Hamiltonian with the extended constraints and the 
extended constraints with each other vanish weakly. The presence of scalar ghosts is 
then exploited to remove the gauge arbitrariness by imposing gauge fixing on them. 

We organize the paper as follows: In section 2 we consider a 1-dimensional 
quantum system to illustrate the supersymmetric extension and to construct the 
functional integral representation of the system transition amplitude in the extended 
space. In section 3 a simple mechanical model with a gauge symmetry is considered 
to show how the ghosts can be added to the theory and how gauge fixing of the 
variables describing the original gauge system can be avoided. We emphasize that 
the choice of simple mechanical models to illustrate the procedure is only for the 
convenience of presentation since the generalization is straightforward. In section 4 
the scheme is applied to Yang-Mills theories. The implementation of the scheme to 
evaluate Green's functions is given in section 5. Section 6 discusses the relation to 
normal gauge fixing and section 7 contains our conclusions. 



2 Ghost extension of a quantum system. 

Consider a 1-dimensional quantum system with Hamiltonian 

Hs = j + V{x), [x,p]=t. (2.1) 

We denote by \s) (or \ip)s) vectors in the system Hilbert space. 
Consider the Hamiltonian 

Hgh = pIpv + pIPz + (^"^{x) {pz + 57^57) . (2.2) 

Here (17^)^ = ff = {pl^Y = = 0, i.e., they are Grassmann canonical operators, 
while z, Pz and their adjoints have boson statistics: 

[^'Pi] = [z\Pz] = i , = -[fl\p^] = i ■ (2.3) 

Here [ , ] is the supercommutator [v4, 5] = AB — {—ly^'^'^ BA where eA,B are 
Grassmann parities of A and B, i.e., e = for bosons and even elements of the 
Grassmann algebra and e = 1 otherwise. We shall call these variables ghosts. 

We denote by \gh)^ (or \i')gfi) vectors in the (indefinite) Hilbert space of system 
( |2.2| ). Note that Hgh depends on a real parameter x, therefore the index x on the 
states. 

We construct a coordinate (Schrodinger) representation of the algebra ( p.3| ) on 
the space of functions ipghii^) = {i^\4')gh = {^y ^* yVyV*\4')gh "with the inner product 

'gkmrgk = j dKC^l^{K)ijgh{^), (2.4) 
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where dn = dz*dzdri*dri. It is easily verified that 

fj = 1] , fj^ = v\ Pn 

Z = Z , = z* , 

in this representation. Note that the dagger denotes the adjoint with respect to the 
inner product ( |2.4[ ) so that the Hamihonian ( |2.2D is self-adjoint. The inner product 
( |2.4| ) is indefinite and the space contains zero and negative norm states. We remark 
that negative norm states always occur when quantizing Grassmann Lagrangians 
bilinear in generalized velocities, which is the case in our fermionic ghost sector (see 
(|2.28| )) as well as in the conventional BRST-BFV formalism. 

The ghost Hilbert space can also be constructed as a Fock space. Define annihi- 
lation operators 

bi = {ujz + ip,)/V2[J , b2 = (coz^ + ipl)/V2^ , (2.7) 
ci = {ujfj + ip^) / V2uj , C2 = {oofj'^ + ipli) / V2uj , (2.8) 

and creation operators through their adjoints. From ( |2.3| ) follows 

[k b\] = [b2, bl] = [ci, c\] = 1 , [C2, cl] = -1 . (2.9) 

Defining the vacuum, |0)g^, as the state annihilated by all annihilation operators 

and with gh{0\0)gfj^ = 1, one notes from (|2.9| ) that states containing cl have negative 
norm. In Fock space the ghost Hamiltonian ( p.2|) is given by 

Hgh = ujNgh = uj(b\bi + b\b2 + c\ci - c\c2) , (2.10) 

where Ng^ is the ghost number operator; it is non-negative, Ngh > 0. 
Define the fermion operators 

Q = c{k - b\c2 , R = c\h + b\c2 (2.11) 

and their adjoints Q'^ and w . We have 

Ngu = [Q,R^] = [Q\R] (2.12) 

and 

[iVg,, Q] = [iV,,, gt] = [Ngr,, R] = [iV,,, R^] = Q . (2.13) 

The ghost Hamiltonian exhibits an iV = 2 supersymmetry generated by Q, R and 
their adjoints. 

Introduce a system with the extended Hamiltonian 

H = Hs + Hgh, (2.14) 
4 
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where Hgh was defined in ( |2.2| ) and u'^ = uj'^{x). Therefore states in the ghost Fock 
space depend on the operator x and accordingly act as operators (functions of x) on 
the system Hilbert space. This motivates the following definition of the extended 
space: 

m = \9hr-\s), (2.15) 

where the dot stands to emphasize that \gh)^ is an operator acting on the system 
state \s). 

In the coordinate representation \k) ■ \x), k\K) = k\k) we have 

{x\ ■ {Km = {x\ ■ {K\ghr ■ \s) = {K\ghr{x\s) (2.16) 

since x is diagonal. 
We observe that 

[H,Q] = [H,Q^] = (2.17) 

since [x, Q] = [p, Q] = 0, but [H, R] ^ as [p, R] ^ (while [x, R] = 0). The N = 2 
supersymmetry is therefore explicitly broken to iV = 1. In the Appendix we discuss 
an = 2 supersymmetric extension of the system dynamics. 

Consider the subspace in the extended space defined by the conditions 

Q\ph)) = Q^ph)) = . (2.18) 

In the Fock space representation of the ghost Hilbert space, one can show that 
the solutions of ( p.l8| ) are 



\ph))=\\0rg,-\s) + J2\<Pnr,h-\Sn). (2.19) 

71=1 

Here |s) and |s„,) are system states and 

I0n)^. = f-io)^, = gl-f-^io)^, = gt<l'f"-i|o),\ (2.20) 

with 

f = bin + c\c\ , $ = h\cl , <!>' = blc\ . (2.21) 

Note that {{ph\ph)) = if A = in (|2.19| ) since vectors (|2.2CI| ) have, by definition 
( |2.18| ), zero norms. 

We remark that the choice of Fock space to discuss the properties of the states 
(|2.18| ) is only a matter of convenience. Indeed, it is easy to prove that these prop- 
erties hold generally, i.e., a state satisfying ( p.l8|) in an arbitrary Hilbert space is 
the sum of a closed state (belonging to kern (Q) /image (Q) n kern((5^)/image((5^)) 
and an exact state (expressible as Q or on some state). By definition ( |2.18| ) the 
latter have zero norm. In Fock space the ghost vacuum is the only closed state. 

For any operator O commuting with Q and it follows from ( ^.2U[ ), ( |2.18| ) and 
(|2l9| ) that 

((p/^'|6b/.)) = (.'|-^,(0|6|0)^,-|.) (2.22) 



if A = exp(ia) and a is independent of x (Even in the case where a is x-dependent 
this relation still holds, provided that p is replaced by p — da/dx). For = 1 this 
yields {{ph'\ph)) = {s'\s). 

If O = Og = 0{x,p) equation ( |2.22D becomes 



{{ph'lOM)) = {s'\Os\s) . 



(2.23) 



To prove ( 2.23| ) we note that since x commutes with 10)^^^, it acts directly on the 
system state. We only have to prove that the same holds for p. From the relations 
(|27^ ), (|2.9|) one easily verifies [p,bi] = phl, where (3 = —iuj'/{2uj), and similar 
relations for the other creation and annihilation operators. From the definition 
of the ghost vacuum we have [p, foilO)^;^] = bi[p, \0)gfj] + /9&2|0)g/i = analogues 
relations for the other annihilation operators. As p commutes with Q and these 
relations imply that the ghost state [p, [0)^;^] belongs the physical subspace ( [2.2U| ) 
with particle number Ngh = 2. We find from (ICTj) [p, |0)^;,] = 0^0)^^,^ = Q\ghy. 
Thus for Os = we deduce 



(2.24) 



Relations (|]2|) as well as (|2l0| ), (l27l^) and {^J§ imply that {{ph'\H\ph)) = 
{{ph'\Hs\ph)) = {s'\Hs\s)- Since both H and Hs commute with Q and their 
application to a physical state leads to a physical state and a simple inductive 
argument then shows that {{ph'\H"-\ph)) = {{ph'\H2\ph)) = {s'\H'^\s) where ( ^.23| ) 
was used in the latter equality. We conclude 



{{ph'\e~''''\ph)) = (s'|e-^*^^|s). 



(2.25) 



As an application of ( p.25| ) we take \s) = \x) and \s') = \x'). This yields 



{{ph'\e-''"\ph)) = I [dx]e 



(2.26) 



where Ss{x) = /g dt^x"^ /2 — V{x)) is the system action. The trajectories contributing 
in the functional integral obey the boundary conditions x(0) = x and x{t) = x'. 

The amplitude in ( |2.26D can be represented as a functional integral over the 
extended configuration space. Consider the functional integral representation of the 
transition amplitude in the extended configuration space 



{{x\K'\e-''"\x,K)) = I [dx][dz*][dz][dri*][dri]e 



(2.27) 



where \x,k,)) = \k) ■ \x) and 



Se = Ss+ / dt'{z*Z + T]*T] - UJ\X){Z*Z + 7]*!^)) =Ss + Sgh . (2.28) 



6 



The trajectories for the ghost variables obey the boundary conditions k(0) = k and 

The functional integral ( p.27|) coincides with ( p.26|) if we calculate its convolu- 
tion with any two physical states ( |2.1(j| ). There is a particular choice of bound- 
ary conditions for the ghost fields in ( |2.27| ) such that ( p.27|) coincides with ( |2.26|) . 
This is achieved by choosing = k = 0. This implies that we choose in ( ^.261 ) 
\ph)) = |0) ■ |x) and \ph')) = |0) • \x') where k\0) = (not to be confused with the 
ghost vacuum). To verify that these states satisfy ( |2.18D , one expresses the operators 
Q and through relations ( |2.8| ), ( p.9[ ) in terms of the ghost canonical variables. 
Using the coordinate representation (|2.5|), (p.6|) we obtain 



{k\Q\0) = {k\Q^\0) =0 (2.29) 

since {k\0) = 5{k) (/ c/K5(K)^/'(/t) = ^'(0)). 

Note that ( ^.27| ) is a Gaussian integral for the ghost variables. Therefore 

J [dn] e'^^"- = F{x, k\ k)Af/Ab = F(x, k', k) (2.30) 

where Ap = As = det^—d^—uj'^) and F{x, k', k) = exp{iSgh{nd)) with a classical 
solution obeying the boundary conditions k(0) = k, and K,{t) = k,'. 
Relation (|2.26|) implies 

J dndn' l^^{ph'\K')F{x, k\ K){K\phYg^ = 1 . (2.31) 

For vanishing ghost boundary conditions, {K\ph)^f^ = 6{k,) and F[x, 0, 0] = 1. 



3 A gauge model. 

In this section we illustrate, by means of a simple model with a SU{2) gauge sym- 
metry, how the method of the previous section can be used to construct a functional 
integral for a gauge theory without fixing a gauge. 

Let the canonical variables p and x be elements of a linear unitary representation 
of SU{2). We choose the system Hamiltonian as in (|2.1| ), but now p'^ = {p,p) where 
( , ) is the invariant inner product. We denote by A" the SU{2) generators in 
the chosen representation. The Hamiltonian (|2.1|) describes a gauge theory if the 
physical states are gauge invariant [|l|, i.e., 

a^lp/i), = ixX''p\ph)s = i{x, X'p)\ph), = (3.1) 

and [HsjCTg] = 0. The canonical commutation relation [x,p] = i implies [(7^,(3"^] = 
2^ga6c^c gQ ^i^o;^ f^i^Q constraints generate SU{2). The Hamiltonian ( |2.1| ) is invariant 
under gauge transformations generated by a", p ^ Tgp, x TgX and TgT^ = 

J-gJ-g -L- 
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As described in section 2, we extend the system Hamiltonian (|2.1| ) by adding a 
ghost Hamiltonian H = Hg + Hgh, where 



Hgh = iPlP.) + {vlVn) + {z\n{x)z) + {f,\n{x)f,) (3.2) 

and f2 is a hermitian, strictly positive matrix. For simplicity we assume the ghost 
variables to realize the fundamental representation of SU{2), i.e., they are complex 
isospinors. 

To make ( p.2|) invariant under the SU{2) transformations z TgZ, fj — > T^r), 
with Tg 9.2x2 unitary matrix in the fundamental representation of SU (2) , we require 
the following transformation rule for Q: Q{Tgx) = Tgfl{x)T^ so that fl transforms 
in the adjoint representation. 

Introducing the fermionic operators 

Q = ^[iv^,Pz) - iPlz)] , Rix) = z[(pt,p,) + iv\nix)z)] (3.3) 

and their adjoints R^{x) we find 

HgH = [Q,R\x)] = [Q\R{x)]. (3.4) 

Note that the operators Q, commute with the total Hamiltonian, while R{x) and 
W{x) commute only with the ghost Hamiltonian. 

As before the Hilbert space of the ghost sector can be build as a Fock space. 
Indeed, if ujf{x) are the eigenvalues of fl{x), then Q = J2iQi and R{x) = J2i ^i{x)Ri 
with Qi and Ri determined by (p.ll|) for each SUSY oscillator labelled by i = 1, 2. 

The constraints are also extended a" = a" + a^f^ with 

agf^ = iz^T^p^ - iplr^z + ifj^^pr, - iplr^fi = [Q, izW^p^ - iplr^fi] . (3.5) 

Here are the Pauli matrices. By construction [a", ^] = 0, [o"", Q] = [a", = 
and a"- generate SU{2). 
We prove the relation 

{{ph'\e-''^\ph)) = s{ph'\e-''^^\ph), , (3.6) 

where the physical states \ph)) and \ph')) in the extended space have non-zero norm 
and satisfy the conditions 

Q\ph)) = Q^ph)) = a^lph)) = 0, (3.7) 



while the system states on the right of ( |3.6D obey ( p.ll) . Relation ( p.6|) establishes 
a one-to-one correspondence between physical transition matrix elements of the 
original gauge theory and its supersymmetric extension. 

From ( |2.25| ) we note that it is sufficient to prove that (|3.7| ) imply ( |3.1| ). Consider 
the second order Casimir operator C = a"" a"'. From ( |3.5| ) it is easily established 
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that C = Cs + [Q, A] where Cg is the second order Casimir of the system, and the 
exphcit form of A can be read of from (|3.5|) . It follows 



{{ph\C\ph)) = {{ph\CM)) = {s\C,\s) , (3.8) 

where ( [^.23D was used. Since the system Hilbert space has a strictly positive norm, 
(|3r^ ) implies (|3.1|) . Thus relation (|3l6|) is proved. 

Next we turn to the functional integral representation of (|3.6|) . To avoid diver- 
gencies in the measure of the functional integral caused by the gauge invariance, we 
have to solve the constraints a"'\ph)) = explicitly by choosing a parameterization 
(coordinates) in the physical configuration space, i.e., we should fix a gauge in the 
extended theory. In contrast with the BRST-BFV scheme we have boson ghosts 
so that we can impose the gauge condition on them, while the system degrees of 
freedom remain free of any gauge condition. 

To fulfil this program, we choose the coordinate representation ( p^.l6| ) {x\ ■ {n\ph)) 
= ipph{x, k), where all states are functions of x and k. To project the total Hamil- 
tonian on the subspace of gauge invariant functions ipph{x^n), we change variables 
in the superspace 

Z = fgXp/\/2, X-^TgX, IJ^fgT], (3.9) 

where = (l^O); P is real and Tg is the group element Tg in the representation 
of X. Denoting the group parameters by 6^, the constraint operators in the new 
coordinates assume the form a" ~ d/dO"' Therefore physical states (gauge 
invariant states) are functions independent of 0". We transform H to the curvilinear 
supercoordinates ( p.9| ) and drop all terms containing d/dO"' to obtain the physical 
Hamiltonian 

Hp,, = Hs + ^ + ^+pIp, + ^{^"f + jx^nx + v^nv , (s.io) 

where a" = a" + aj with aj = ivj^r'^p^ — ipl^T""!] and Pp = —ip^^^'^dp o p^/^ is 
the hermitian momentum conjugated to p. The third term in ( 3.10|) is due to the 



operator ordering. Restoring Planck's constant, this would be proportional to h^. 

The above procedure to obtain the physical Hamiltonian also applies to the 
operators Q and to derive the supersymmetry generators on the space of functions 
invariant under the gauge transformations generated by a". We find 

X z 1 

Ql = --^xVp + ^x^p.-^[xV + ^W + ^^')] (3.11) 

with ip"^ = (0, 1). Taking into account that t] and r/^ do not commute with a'*, one 
can verify that these two operators are adjoints with respect to the measure p^dp. 
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Indeed, the term proportional to 1/p, which appears when the adjoint of dp is taken 
with respect to this measure, precisely cancels against a similar term originating 
from the reordering of r] or t]'^ and a'^ in the adjoint of the third term. It can also 
be checked explicitly that these operators commute with the physical Hamiltonian 
(|3.1CI| ) as they should. 

The physical states of the system are obtained as before by imposing 

QM)) = Ql\ph)) = Q. (3.12) 

We note that the system variables x and p do not commute with the supersym- 
metry generators Qp and Q],- Indeed, the operator S}Qp + Q\i, with and ^'^ 
Grassmann variables, involves a linear combination of the system constraint opera- 
tors a^. Therefore it generates an isotopic rotation of x and p of the same form as 
an infinitesimal gauge transformation, but with the parameters being even elements 
of a Grassmann algebra. Therefore the condition ( |3.12| ) ensures gauge invariance in 
the vacuum sector of the ghosts, i.e., a vector describing any system excitation will 



automatically satisfy the Dirac condition (p.l|) . 



Note how the reduction in the degrees of freedom is obtained: we start with 
n system degrees of freedom, n being the real dimension of the representation in 
which the x transforms. Then we add 4 real bosonic and 4 real fermionic ghost 
degrees of freedom in a supersjnnmetric combination. Since the system variables 
realize a trivial representation of the supersymmetry an exact cancellation of the 
ghost degrees of freedom is ensured, leaving only the n system degrees of freedom. 
Next we eliminate three real bosonic ghost degrees of freedom by projecting on 
the subspace of gauge invariant functions (see ( p.9|) ). The projection (|3.9| ) mixes 
the system and ghost degrees of freedom in a non-linear way. In particular the 
angular variables of the bosonic ghosts are absorbed in the new system variables so 
that three angular variables of x are the supersymmetric partners of three fermionic 
ghost degrees of freedom. The condition ( |3.12| ) ensures that the 4 fermionic ghost 



degrees of freedom are balanced by the p ghost degree of freedom and 3 system 
degrees of freedom. Hence, we are left with — 3 degrees of freedom, which is 
correct as we have 3 independent constraints ( p.l|) in the original model. 

In our arguments above we have assumed that the system constraints a" are 
irreducible, i.e., they are all independent, which would not be the case if, for example, 
we take x to be in the adjoint representation. Nonetheless, in our approach the 
supersymmetry condition (|3.12|) will still provide the elimination of the right number 
of unphysical degrees of freedom in the system sector, provided the ghost sector is 
chosen so that the constraints a^f^ ( |3r5| ) are irreducible. The latter is always possible 
so that, in contrast to the conventional BRST-BFV formalism, we do not need ghosts 
for ghosts if the system exhibits reducible constraints. 

Returning to the functional integral representation, let Hp denote a set (p, r/^, rj). 
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Using the formalism of [Q we find 



((x', K'\e-''^^-\x, Kp)) = [ dK''dx" ip'p"')-'/W^ix', k' x", k'')Qs{x'\ k'' X, 



(3.13) 

the amphtude Ui^^ is given by a functional integral for the extended theory in the 
unitary gauge z = XP- 



Here 



and 



U:^fix,K,,x",K;) = j[dp][dx][dK^][dp,^]e''^» (3.14) 

- J[dx][dy][dK,][UUo\p\']e''''\ (3.15) 

= I dt'{xp + pI^t] + r]^pr^ + ppp - Hph) (3.16) 

'J 

1,„ ,^ . 1 



-Vix) - -p\^nx - V^nr] - — ] . (3.17) 



We denote by Dt = dt + iy°'\°' and Dt the covariant derivatives in the representation 
of X and the fundamental representation (Aa Tq), respectively. 

To go from ( p.l4|) to ( p.l5| ) we first replace the measure [rf^p] [rfp^p] by the 



Faddeev-Popov measure in the unitary gauge [dK\[dpi^Y\.t^FP^{z — xp)'5(o""). The 
Faddeev- Popov determinant for the unitary gauge z = px is App ~ \p\^. Corre- 
spondingly in (|3.14|) is replaced by the Hamiltonian action for the total extended 
Hamiltonian. Next the identity J[dy°']exp{iy°-a°') = Ilt6{a°') is used and the inte- 
grals over all momenta are done. Finally the integral over z and are done using 
the delta function and (|3.15| ) is obtained. Note that apart from the operator order- 



ing term (~ h / p^) in (|3.17|) all the other terms depending on p can be written as 



DtA'^ ^ z^^z where z = px- Thus ( p.l7| ) is the action of the extended supersym- 



metric gauge theory in the unitary gauge imposed on the boson ghost z. 

The kernel in p.l3| ) takes care of the global properties of the change of 



variables H, P] 

Qs{x",KpX',K'p) = "^Si^x" - Sx')6gh{Kp, SKp) (3.18) 
s 

(/ dn'Sghi^', = V^(^))- The sum is over all inequivalent elements s acting on 

X and Kp as follows: sx = TgX and SKp = TgUp where, by definition, TgX = ±X ci-nd Tg 
are corresponding SU{2) elements in the representation of x. The s-transformations 
form a Z2 subgroup of SU{2), Tg = ±1. 

We note that any physical state ip{xp, Kp) is invariant under the s-transforma- 
tions. Here we denote the transformed system variables defined in ( p.9|) by Xp 
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to distinguish from the original system variables. We also define the action of s 
on 6*" by Tg{s6) = Tg{0)T~^. Then sx = x and = k so that ip^sx^sn) = 
ip{x, k) for any element of the extended Hilbert space. For a physical state we have 
ipph{x,K) = ipph(Tg{9)xp,Tg{9)Kp) = ipph^Xp, Kp). Since s leaves x and k invariant, 
we can replace 6°", Xp and Hp by their s-transforms in the latter equality. Thus 
we conclude ipph{sxp, SKp) = ipph{xp, Kp). The time evolution must respect this 
symmetry. However, the functional integrals (|3.14|) or ( p.l5|) do not respect this 



symmetry so that they cannot describe the correct time evolution. The particular 
linear combination of functional integrals ( |3.15| ) with different boundary conditions 



as prescribed by (|3.13|) and ( |3.18|) has the desired invariance property @, |T3| 

From ( |3.6| ) we know that ( |3.13| ) coincides with the physical system transition 
amplitude only after supersymmetric and gauge invariant boundary conditions have 
been imposed on the ghosts. In analogy with section 2, the state {{x,K\ip)) = 
Sgh{K)ipsix) is supersymmetric, i.e., it is annihilated by both supersymmetric gen- 
erators. However, it is not gauge invariant so that it cannot be used as boundary 
conditions in ( p.l3| ). To construct a gauge invariant supersymmetric state we average 
5gh{n)ips{x) over the gauge transformations. This yields 

{{Kp,x\ph)) = ^S{p){v^vrrs'ix) = -^S{z^z){v^vfrs\x) , (3.19) 

where ip^^lx) is any gauge invariant function of x. The ghost part of ( p.l9| ) is 
just a delta function with respect to the scalar product measure dpp^ J dri'^drj. 
Substituting ( |3.19| ) into ( p. 6] ) and using ( |3.13| ) we obtain 

s{ph'\e-''^'\ph), = J dxdx'{ph'\x')Uf\x',x){x\ph) (3.20) 

with 

U!"ix\x) = hm (ppT'^W^ix'^K'x^Kp) . (3.21) 



There is no problem to impose zero boundary conditions on the fermion ghosts in 
the functional integral ( p.l5| ) . For the boson ghost p one can only prove that the limit 
( 3.21| ) exists. Indeed, in the neighborhood of p = or p' = 0, the integral ( glSl ) 



satisfies the Schrodinger equation {—dp/2 + 3/(8p^))?7j = + 0(1). Therefore 
U^^^ ~ (pp')'^/^ as p and p' approach zero and the limit ( p.21| ) exists. 



4 Yang-Mills theories. 

Here we briefiy outline how the procedure of the previous sections can be generalized 
to Yang-Mills theories. A more detailed discussion can be found in |l^ . 



One can go from the gauge model of section 3 to the Yang-Mills theory (with 
SU (2) gauge symmetry) by simple analogy. We identify the vector gauge potentials 
Ai , (i = 1, 2, 3) with the system variables x, the Lagrange multiplier y with Aq and 
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K becomes a set of ghost fields. One can convince oneself that the expressions ( p. 201 ) 
and ( |3.21| ) still hold with the replacements x ^ Ai and dx 'n.xdAi{x), where the 
product is taken over all space points. Similarly the measure is replaced by the 
functional measure /i = IlxP^ix) and 

Uff[A\ A, = j [dA,][dp]Mdri^][dri] e^''" (4.1) 

with 

Here is the covariant derivative in the fundamental representation and a specific 



choice of the matrix VL (see ( |3.2|) ) was made, Vt = —Df{A) + m^, (in order for Q to 
be strictly positive a mass term is added for the ghosts). The quantum corrections 
to the potential energy is proportional to a divergent factor Vq{p) = 3[(5^(0)]^/(8p^) 
as is usual for operator ordering terms in a field theory. Making use of a lattice 
regularization of space in ( [4.2|) our earlier arguments can be repeated to prove the 
existence of the limit ( |3.21| ) for each degree of freedom of the field p{x). The coupling 
of the field oscillators due to the term (9jp)^ is of no consequence in the limit p — 
as it is negligible in comparison with the singular quantum potential Vq. Therefore 
the limit ( 3.21| ) exists and supersymmetric gauge invariant boundary conditions can 



be imposed as before. 

Due to the presence of the singular potential Vq{p) the limit (p.21| ) is highly un- 



desirable and leads to technical difficulties. This limit originates from our choice of 
supersymmetric boundary conditions for the path integral. Fortunately, all physi- 
cally relevant information about a quantum system can always be extracted from 
Green's functions. In the next section we show that the restriction of the supersym- 
metric boundary conditions can be dropped in the path integral representation of 
Green's function and, hence, the aforementioned technical difficulty can be avoided. 



5 System Green's functions. 

The quantities one is normally interested in are the Green's functions of the system 
from which the S-matrix elements and the excitation spectrum can be extracted. 
They are defined as the vacuum (groundstate) expectation values of time ordered 
products of the field operators. For complex time they are defined by the following 
quantities, taken in the limit r/ — * oo, Tj — > — oo: 

{s'\e-^^^fdiin)d^,iT,) . ..d^^iT^)e^^%s) 

Here r/ > ri > . . . > r^, \s), \s') are arbitrary system states and 0|(ri) = 
^HsTiQs^-HsTi g^j^g system operators in the Heisenberg picture. The real time Green's 
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functions are obtained from the complex time Green's functions by analytic contin- 
uation (Wick rotation). Following the standard procedure one can also express the 
Green's functions as functional integrals. 

When the system possesses a gauge symmetry, it is again necessary to eliminate 
the unphysical degrees of freedom to obtain well defined Green's functions. Follow- 
ing the Faddeev-Popov approach, one is again faced by the difficulties associated 
with gauge fixing. Therefore we would like to extend our procedure to the evalua- 
tion of system Green's functions. In principle this is straightforward, except for the 
following point which require special care: The groundstate of the supersymmetric 
extended Hamiltonian may turn out not to be supersymmetric, i.e., the supersym- 



metry is spontaneously broken and the groundstate does not satisfy ( |2.18| ) and is 
thus unphysical. If this happens the Green's functions of the supersymmetric ex- 
tended system do not coincide with those of the original system, as the expectation 
value of the system operators are evaluated with respect to an unphysical state. One 
way of avoiding this difficulty is to impose supersymmetric boundary conditions in 
the functional integral for the Green's functions as relation (p.25|) holds regardless 
whether the supersymmetry is spontaneously broken or not. From our previous dis- 
cussion it is, however, clear that these boundary conditions are inconvenient from 
a calculational point of view. Instead we present here an alternative way of ob- 
taining the system Green's function which allows us to choose arbitrary boundary 
conditions. 

For the purpose of illustration let us consider again the 1-dimensional quantum 
system of section 2. The complex time Green's functions of the supersymmetric 
extended system is obtained from a simple generalization of ( p. ID 

((V''|e-^(-/--»)|V)) 

Here \ip)) = \gh) ■ |s), l^j')) = \gh') ■ \s') are states in the extended Hilbert space and 
H is the supersymmetric extended Hamiltonian 

H = Hs + aHgh (5.3) 

where Hs, Hgh were defined in ( p.l| , p.2|) and a is an arbitrary real number. The op- 
erators Oj(rj) = e^'^'Oje"^'^' are supersymmetric extensions of the system operators 
in the Heisenberg picture defined by Oj = + [Q, (Oj)^/,] with any desired choice 
of {di)gh. 

Note that the extended Hamiltonian above differs slightly from that of (p.l4|) in 



the presence of the arbitrary real number a. Since the ghost Hamiltonian is irrelevant 
on the physical subspace (see (|2.23| , |2.25|) ), physical quantities will be independent 



of a, while excitations in the ghost sector will be affected. We now exploit this 
fact to our advantage by considering the limit a ^ 00. From eq. ( |2.10|) we note 
that this limit will push all states with non-zero ghost number to infinite energy, 
while the states with zero ghost number will be untouched. Thus in this limit the 
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groundstate of the extended Hamiltonian will be supersymmetric and assumes the 
form \0)gf^ ■ \gs)s, where \gs)s is the system groundstate (see ( |2.19| )), so that the 
Green's functions of the extended Hamiltonian coincide with those of the system 
Hamiltonian. We now prove that this is indeed the case. 

Let \N,x)) = \N)^f^ ■ \x) be eigenstates of Ngh and x where we noted that 
[Ngh,x] = 0. The states \N,x)) form a complete basis in the extended Hilbert 
space {\N)gf^ is the Fock basis of the ghosts). Writing out the Heisenberg operators 
in (|5.2|) and inserting a complete set of states in the form 

l = J2 f dx\N,x)){{x,N\ (5.4) 

before and after each time evolution operator, the Green's function ( ^.2| ) reduces to 
the product of generic matrix elements of the extended system operators ((x', A^', |0j 
\N,x)) and the time evolution operator {{x' ,N'\e^'^^^^^"-^ah)\]\j^x)) = {{x\N'\Ut 
\N,x)) where T = tj — Tj+i > 0. Consider the latter matrix elements. From 

Hgh\N,x)) = uj{x)N\N,x)) (5.5) 

follows that in the limit a — * oo a perturbation expansion in Hg can be developed, 
using aHgh as the unperturbed Hamiltonian. Introducing an interaction picture we 
can express the time evolution operator as: 

oo T Tl r„_i 

Ut = e-"^«'^^^(-l)"y'rfriy'rfr2... 1 cir„^,(ri)^,(r2) . . . i^,(r„) 



n=0 







E(-ir^n, (5.6) 



where T > ti > r2 . . . r„ > and i/s(rj) = e^^ahnfj^f^-aHghn^ 'Jq evaluate the 
matrix elements {{x' , N'\Ut\N, x)) , consider a typical term in the series ( |5.6|) and 
insert the complete set of states (|5.4| ) between the Hs{Ti). Using ( ^.5|) we find: 

{{x\N'\B^\N,x)) 

T n -r„_i 

= j dxi...dxn-i JdnJdT2... J cir„e-"'^'^'(^-^i)((x',iV'|i^,|iVi,xi)) 

Ni...N„-i 

xe-'^'^'^'^^'-^^^{xi,Ni\Hs\N2,X2)) . . . g-^'^-i^-i^^-i-^-) 

X {{xn-1, N^-i\Hs\N, x))e-"-^^" . (5.7) 

Here the shorthand notation u' = uj{x'), uJi = uj{xi) and u = u!{x) was used. 
Since T > ri > r2 . . . r„ > we deduce that in the limit a oo only terms with 
N' = Ni = . . . = N = survive. The ghost vacuum (A^ = 0) is a physical state 
and therefore ((x', 0|ifs|0, x)) = {x'\Hs\x) according to (|2.23|) . Making use of these 
observations we find that eq. (|5.7|) becomes 



{{x',N'\B^\N,x)) ^ —{x'\H:\x)Sn',oSn,o, (5.8) 
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when a ^ oo. From (|578|) we conclude 

{{x',N'\Ut\N,x)) {x'\e-^'^\x)5N',oSN,o, a ^ oo . (5.9) 
Substituting this result in (|5.2| ) and using ( p.23|) , we obtain in the limit a — oo 



This establishes that by extending the Hamiltonian as in ( |5.3| ), and taking the 
limit a oo before the limits tj —>■ oo, Xj — oo, we indeed recover the system 
Green's functions from the Green's functions of the supersymmetric extended sys- 
tem. In this approach arbitrary boundary conditions can be imposed in the path 
integral for the Green's function as the limit a — > oo ensures that the supersymmet- 
ric groundstate is selected in the evaluation of the expectation value. 

We conclude this section with a couple of general remarks. The first remark 
concerns degeneracies. We note that for some values of x, x', N and A^' the following 
equality holds: uj{x)N = uj{x')N' and we have to analyze what happens to the 
perturbation expansion ( |5.7| ) when this degeneracy occurs in the ghost energy levels 
(see (|0| )). This degeneracy does not affect our conclusions (^.8|) and, hence, (|5.1CI|) . 



Indeed, the operator is local and therefore its matrix elements {x'\Hs\x) have 
support only on x = x'. Consequently this equality only holds when N = N' and it 
is necessary to consider the degeneracy of a ghost level with a fixed ghost number 
only. In this case some of the exponential factors disappear and the corresponding 
integration yields a typical factor TigN, with the degeneracy of the state \N)gh. 
Since this is a finite factor, independent of a, we conclude that in the limit a oo 
the factor e""'^^^ ensures that only contributions from the ghost vacuum survive. 

The above argument will hold for any supersymmetric extension of the Hamil- 
tonian, provided that the ghost Hamiltonian itself does not exhibit spontaneous 
breaking of the supersymmetry. 

Clearly the above procedure would not be necessary if one can prove that the 
supersymmetry is not spontaneously broken for a = 1. If the groundstate of the 
extended Hamiltonian has the form ( p.l9| ), it is not hard to see that the Green's 



functions ( [5. 21 ) coincide with the system Green's functions ( |5.1D in the limit Tf ^ oo, 
Ti —OO. However, a proof that the supersymmetry is not spontaneously broken 
does not seem easy. 

We remark that the procedure as outlined above does not apply to relativistic 
systems as it would break Lorentz invariance. However, a relativistic invariant 
formulation of the above procedure can be obtained by taking the limit where the 
mass of the ghost fields go to infinity, m — oo (see ( [4.2|) ). In this limit contributions 
of ghost excitations to the Green's functions are exponentially suppressed because 
their energies grow as m. It is then not difficult to generalize the above argument 
to show that our conclusion (|5.1CI|) holds in this case. 
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In the case of a gauge system, one is interested in the Green's functions of gauge 
invariant system operators. They are obtained from our prescription above by con- 
sidering only physical operators, Oi(rj), in (|5l^) , i.e., operators that commute with 
the supersymmetry generators obtained after projecting onto the space of gauge 
invariant functions. For the gauge model of section 3 we would therefore consider 
operators that commute with the supersymmetry generators Qp and Qj^ of ( |3.11| ). 
From the discussion below ( p.l2|) it follows that any system operator which com- 
mutes with Qp and Qp has to commute with the system constraint operators a", 
and is thus gauge invariant. Therefore the Green's functions obtained from our 
description above coincide with the system Green's functions of gauge invariant 
observables. 

Finally we remark that the evaluation of the Green's functions in a gauge theory 
using the supersymmetric quantization is not an easy task. Though the problem of 
the supersymmetric boundary conditions has been avoided, which is a great simpli- 
fication from the calculational point of view, there is still a lack of a well elaborated 
techniques in quantum field theories to treat, even in perturbation theory, local 
measures in the functional integral (see (|3.15| )) and singular quantum potentials 
(see ( |3.17|) ) that usually occur in unitary gauges. Nonetheless it is possible to cal- 
culate the renormalized Green's function using a modified loop expansion [p 



6 The relation to gauge fixing. 

Here we explain how the functional integral ( [4.1D is related to the functional integral 
of a gauge fixed Yang-Mills theory. For simplicity we consider the model of section 
3. Let X = xp he a. generic configuration satisfying the gauge condition F{x) = 0. 
We assume that the gauge condition fixes all continuous gauge arbitrariness, i.e., 
the stationary group of x = Xi;' is trivial. However, the equation F{x) = may have 
many solutions related by discrete gauge transformations (Gribov copies). This 
gauge arbitrariness does not decrease the number of physical degrees of freedom, 
but it does reduce the volume of the physical configuration space spanned by xp- 
If all configurations xp form an Euclidean space iR*^ and S* is a set of residual 
discrete gauge transformations, then the physical configuration space (a fundamental 
modular domain) is isomorphic to Aj? ~ / S C . 

Consider the change of variables Xp = TgXp, z = TgXP, rj = Tgrjp, where the first 
relation determines the element Tg for a given xp and Xp, and Tg is this element in 
the fundamental representation. The latter two relations define the new variables 
z and f]. Introducing this change of variables in the Hamiltonian ( |3.10| ), the ghost 
part assumes the form ( p.2|) , while in the system Hamiltonian the Laplace opera- 
tor p"^ must be written in the curvilinear coordinates Xp = TgXp. If there were no 
ghosts in the theory, the Laplace-Beltrami operator p"^ would contain operators a^, 
corresponding to momenta conjugate to the group parameters 6"". Since the above 
change of variables involves the ghost variables, all operators a" get replaced by 
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dgf^ just as the boson ghost constraints were replaced by a"^ + in ( CT )- The 
functional integral for this Hamiltonian takes the form ( p.l3| ) where the measure 
is replaced by fi^xp) {J^ dhipp^ J dxp = J dn J^^ dxFp{xF))- The measure pi^xp) is 
the Faddeev-Popov determinant in the gauge chosen. The kernel ( |3.18| ) contains a 
sum over all elements Tg satisfying F{TsXp) = 0. Imposing supersymmetric bound- 
ary conditions for the ghosts by choosing {{K,XF\4'))ph = ^ ghi^^i^f^ {x p) , the ghost 
integral is Gaussian and can be done explicitly. Boson and fermion determinants 
cancel and the prefactor is 1 due to the supersymmetric boundary conditions (see 
( |2.31| )). We therefore end up with the functional integral for the system action in 
the gauge F{x) = 0. 

For the Yang-Mills theory we can choose F{A) = diAi = (Coulomb gauge). 
The above procedure yields the functional integral obtained in p. 



7 Conclusions. 

The normal procedure of Lorentz covariant gauge fixing has the disadvantage that 
it leads to Gribov copying. The description of these copies and the fundamental 
modular domain is extremely difficult and in the continuum theory not well founded 
since it depends strongly on the functional space chosen for the gauge potentials 
(the normal choice of forms a set of zero measure in the functional integral) 
However, to calculate the non-perturbative Green's functions in covariant gauges 



one needs all copies for a generic configuration satisfying the gauge condition 
Given the complexity of these copies the non-perturbative evaluation of Green's 
functions in covariant gauges seems to be a hopeless task. 

In the procedure developed above we have solved these difficulties. Our effective 
theory is Lorentz invariant and at the same time the Gribov problem has been 
avoided. The price we pay is in the appearance of an additional scalar ghost field. 
Furthermore we have encountered a difficulty with the supersymmetric boundary 
conditions to be imposed on the functional integral. This difficulty can, however, 
be avoided in the path integral for the Green's functions as we have shown. 

We would like to stress that in the asymptotic free domain where perturbation 
theory is valid, the normal approach of Lorentz covariant gauge fixing ||^ is more 
appropriate. However, in the infra-red limit where the non-perturbative aspects 
of the theory begin to dominate, the present formalism should come into its own 
right. One of its possible applications would therefore be to study the infra-red 
behavior of the theory. As the present approach is gauge independent it makes it 
possible to extract the infra-red properties of the theory due to the true dynamics of 
the system (self-interactions of the gluons) and free from any infiuences coming from 
the presence of the gauge dependent Gribov horizon in a gauge fixed approach. This 
would hopefully provide insight into the physical meaning of the (non-perturbative) 
dynamically generated mass scale found in 0, |TD|, |TB[. 
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Appendix. 

There exists an iV = 2 super symmetric extension of the system dynamics. To obtain 
it, one should find a ghost extension P — oi the momentum p such that 

[P, R\ = [P, W\ = [P, Q\ = [P, Qt] ^ (A.l) 

and replace the p by it in all system operators. One can convince oneself that the 
operator 

P^pt^p- ip{x){c2C, - clcl + hh - b\bl) , (A.2) 

with /3 = — a;'/(2a;), is the desired N — 2 supersymmetric extension of the momen- 
tum operator. An advantage of this extension is that we know the structure of the 
eigenstates of the extended Hamiltonian 

\E)) = \N)t, ■ \E^)s , (A.3) 
where \N)gf^ are the eigenstates of the ghost number operator Ng^ and 

{H, + uj{x)N)\E^)=E^\E^),. (A.4) 

The N = 2 supersymmetry is not spontaneously broken in this the 
groundstate \gs)) = \0)gj^ ■ \gs)s, with \gs)s the system groundstate, is manifestly 
N = 2 supersymmetric. 

A disadvantage of this extension is the non-locality of the extended system 
Hamiltonian that arises in the field theory case. Recall that is replaced by 

—D'f{A) + w? and, hence, 13 would involve the non-local operator {—Df{A) + im?)~^ . 
For this reason we did not use this extension to resolve the problem with the super- 
symmetric boundary conditions. 
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